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Abstract 

Single variable hypergeometric functions p F q arise in connection with the 
power series solution of the Schrodinger equation or in the summation of per- 
turbation expansions in quantum mechanics. For these applications, it is of in- 
terest to obtain analytic expressions, and we present the reduction of a number 
of cases of p F p and p +iF p , mainly for p = 2 and p = 3. These and related series 
have additional applications in quantum and statistical physics and chemistry. 
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Introduction 

The types of single variable hypergeometric series considered here appear in nu- 
merous physical contexts including quantum chemistry the development of few- and 
many-body wavefunctions, and statistical physics [6, 13, 14, 8]. A particular in- 
stance in statistical mechanics is a study of phase transitions in the n-vector model 
[8]. In addition, they have arisen, for instance, in condensed matter physics in the 
consideration of the effective vortex mass in continuum models of superconductors 
[4]. Hypergeometric series play an important role in expressing Clebsch-Gordon and 
Racah coefficients (3nj symbols) [13, 14]. We recall that the orthogonal Hermite and 
Laguerre polynomials that appear often in standard quantum mechanical problems 
are instances of the confluent hypergeometric function i F 1 . 

Besides these fairly general physical applications, hypergeometric functions of the 
form q F q or q +iF q are important in quantum mechanics in the series solution of the 
Schrodinger equation [6] and in summing perturbation expansions [16]. We give a 
few details and offer some remarks that may serve to unify earlier treatments. 

In Ref. [16], Saad and Hall summed perturbation series for the spiked harmonic 
oscillator, a linear oscillator with Hamiltonian perturbed by an inverse power-law 
term. In their approach, in order to perform sums with summands containing either 
the function ± Fi or 2 Fi, they needed explicit expressions for the function 3 F 2 to insert 
into an integral representation. Among their key intermediate results is the expression 
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([16], Lemma 2) 

z 3 F 2 (a + 1, 1, 1; 2, 2; z) = ^ £ * ^ dt, a ^ 0, |z| < 1. (1) 

We find that if we perform an integration by parts in this equation and then a change 
of variable, 

z 3 F 2 (a + 1, 1, 1; 2, 2; *) = -z £ ^j^ dy, (2) 

a result that is a special case of a representation found significantly earlier by Gottschalk 
and Maslen [6] in connection with the power series solution of the Schrodinger equa- 
tion. In addition, the case of a = 1/2 in Eq. (2) was known previously [9]. 

On the other hand, it is still possible to extend some of the integration results of 
Gottschalk and Maslen [6] for the function 3 F 2 and we do so in a later section. These 
results have a richness of special cases. 

In Ref. [5], Exton recently investigated reductions of the Kampe de Feriet double 
hypergeometric function [3] and found some new relations for single hypergeometric 
functions. These reductions arise when the two arguments are related by x = ±y and 
are essentially the result of the application of summation formulae for nearly poised 
hypergeometric series. In the following section we describe the particular character 
of the single variable hypergeometric function special cases of Exton, and point out 
connections with earlier work. We provide direct proofs of these special cases, and the 
technique may be useful to other researchers in the variety of physical and engineering 
science that we have mentioned. 
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We then give a more general 3 F 2 reduction in a succeeding section. After com- 
pleting and illustrating two 3 F 2 reductions begun by Gottschalk and Maslen [6], we 
supply brief concluding remarks,. 

Hypergeometric reduction formulae 

Exton [5] deduced the following four relations for single variable hypergeometric 
functions of the form q F q or q +iF q : 

e y 2 F 2 (a, 1 + a/2; a/2, b; -y) = 2 F 2 (2 + a - b, b - a - 1; b, 1 + a - b; y), (3) 

(l-y)- h 3 F 2 [h, a, l+a/2; a/2, b; y/(y-l)] = 3 F 2 (h, 2+a-b, 6-a-l; 6, l+a-6; y), (4) 
2 Fi(a, -e/2; l+a+e/2; x) = (l-x) e ^(a+e, l+a/2+e/2, e/2; l+a+e/2, a/2+e/2; x), 

(5) 

and 

2 Fi(a+d-l, -2a; d; y) = (l-y)~ a 4-Ps(a+d-l, -a, 2-d-2a, d+2a-l; tZ, rf+a, l-d-2a; y). 

(6) 

In order for these functions to be well defined, we require that no denominator pa- 
rameter be zero or a negative integer. An examination of Eqs. (3)-(6) indicates their 
special character. In each case, there is one (or more) hypergeometric function(s) 
with a numerator parameter exceeding a denominator parameter by precisely one, 
and this is exploited in the direct proofs. 

Indeed, equations (3)-(6) are partially based upon the identities [17, 5] 

tp / /o,i /(i 1 i\ (2 + a-6) w (6-a- l) n 

3 F 2 (a,a/2 + l,-n;a/2,6;l) = —— — — , (7) 

{b) n (a -b + l) n 
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3 F 2 (a, 6, -n; a - b + 1, 26 - n + 1; 1) = ) , Jn \' - n , , (8) 



(a-26) n (o/2-6 + l) ra (-6) T 
(a-6 + l) n (a/2-6) n (-26) n ' 



and 



4 F 3 (a, a/2 + 1, 6, -n; a/2, o - 6 + 1, 26 - n + 1; 1) = ( ° 6 f jy.^ - ( 9 ) 

As each of the hypergeometric functions on the left sides of Eqs. (7)-(9) has a 
negative integer numerator parameter, they are individually terminating series at 
unit argument. 

Equations (4) and (5) are proved in Slater's book [17] (Section 2.4.2) by means 
of Vandermonde transformations. As mentioned there, Eq. (5) is essentially due to 
Bailey. 

For the direct proofs of Eqs. (3)-(6), we generally transform each side into a com- 
mon intermediate expression. Such alternative intermediate expressions are likely of 
independent interest themselves. As these transformed expressions have fewer numer- 
ator and denominator parameters, they may be suitable for numerical computations. 

Proof of Eq. (3). For the right side of this equation, we have 

2 F 2 (2 + a-b,b-a-l;b,l + a-b;y) = 1 F 1 {b-a-VMy)-\ iF^b-a; 6+1; y) (10a) 

b 

= e^F^a + 1; 6; -y) - V - ^(a + 1; b + 1; -y), (106) 

where i F 1 is the confluent hypergeometric function. This function results from a 2 F 2 
function when a numerator parameter equals a denominator parameter. In obtaining 
Eq. (10a), we wrote the defining infinite series for 2 F 2 , used the relation (a + l)j = 



(1 + j /a)(a)j, manipulated the series, and applied the relation (a)j+i = a(a + 
Here, as usual, (a) n = F(a + n)/T(a), where T is the Gamma function, denotes 
the Pochhammer symbol [2]. In obtaining Eq. (10b), we applied Kummer's first 
transform 1 F 1 (a; p; z) = e z ±Fi(p — a; p; — z). By again using (a+l)j = (l+j/a)(a)j, 
we may rewrite Eqs. (10) as 
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2 F 2 (2+a-b,b-a-l;b,l+a-b;y) = e y [ 1 F 1 (a;b--y)--y 1 F 1 (a + l;b+l;-y)]. (11) 
Similarly, by using specifically 

(l + o/2),-= 1 + 
in the infinite series form of the left side of Eq. (1), we have 
e y 2 F 2 (a,l + a/2;a/2,b; -y) = e y 



a/2 



(a/2), 



(12) 



= e^Fxfa; b; -y) - -y x F x {a + 1; b + 1; -y)\. (13) 

The equality of Eqs. (11) and (13) shows that Eq. (3) holds. 
Proof of Eq. (4). For the right side of this equation, we have 

3 F 2 (h,2 + a - b, b - a - 1; b, 1 + a - b; y) — 2 Fi(h, b - a - 1; b;y) 

~V2F 1 (h + l,b-a;b + l;y). (14) 
By applying a standard transformation of the Gauss hypergeometric function, 2-^1(01, P; 7 
(1 - z)~ a 2 -Fi[a,7 _ P',T, Z /( Z - 1)]; we nave 



3 F 2 (h,2 + a - b,b - a - l;b,l + a - b;y) = (1 - y)~ h 2 F 1 [h, a + 1; b; y/(y-l)] 



= (i-y) 



-h 



~y(l - y)-^ 2 F 1 [h + l,a + l;b+l; y/(y - 1)] 
2 F 1 [h, a; b; y/(y - 1)] - y 2 F 1 [h + 1, a + 1; 6 + 1; j//(y - 1)] 



(15) 



On the other hand, for the left side of Eq. (4), by using Eq. (12) we have 



(1 - y)~ h sF 2 [h, a, 1 + a/2; a/2, b; y/(y - 1)] = (1 - y)~ h [ 2 F 1 [h, a; 6; y/(y - 1)] 
, 2a y 



2 F 1 [/i + l,a + l;6+l;y/(y-l)] 



(16) 



b (y-1) 

The equality of the expressions in Eqs. (15) and (16) verifies Eq. (4). We note 
explicitly that Exton's Eq. (13) [5] on the left side has omitted a — sign in front of 
the argument y/(l — y). 

Proof of Eq. (5). We proceed similarly to the above proofs, although now we 
require contiguous relations [2, 7, 17] for the function 2 F 1 . We first transform the 
function 3 F 2 appearing on the right side of Eq. (5): 

3 F 2 (a + e, 1 + a/2 + e/2, e/2; 1 + a + e/2, a/2 + e/2; x) = 2 F 1 (a + e, e/2; a + e/2 + 1; x) 



a + e/2 + 1 



x 2 Fi(a + e + 1, e/2 + 1; a + e/2 + 2; x). 



(17) 



We then apply 



,Ma, (3; r , z) = (l- zy- a -P 2 Fx(7 - a, 7 - /3; r ,z), 



(18) 



giving 



3 F 2 (a+e, l+a/2+e/2, e/2; l+a+e/2, a/2+e/2; x) = (l-x)- £ [(1 - x) 2 Fi(a + 1, 1 - e/2; e/2 + 1; x) 



+ 6 x 2 F 1 (a + 1,1- e/2; a + e/2 + 2; x) . (19) 

We use two different contiguous relations (e.g., [7]) to re-express each of the 2-fi's 
on the right side of Eq. (19). We have 

a(l - x) 2 F 1 (a + 1,1- e/2; a + e/2 + 1; x) = (a + e/2) 2 Fi(a, -e/2; a + e/2; x) 

~ 2 F 1 (a, 1 - e/2; a + e/2 + 1; x) 
= (a + e/2) 2 Fi(a, -e/2; a + e/2; x) - | 2^1(0, -e/2; a + e/2 + 1; x) 

x 2 F ± (a + 1,1 -e/2; a + e/2 + 2; x). (20) 



2(a + e/2 + l) 

For the second 2 F\ 011 the right side of Eq. (19) we have 

~ax 2 Fi(a+l, l-e/2; a+e/2+2; x) = (a+e/2)(a+e/2 + l) [ 2 -Fi(a, -e/2; a + e/2; x) 

- 2 Fi(a,-e/2;a + e/2 + l;a:)]. (21) 

This equation is also used to replace the last term of Eq. (20) with x 2-^i(---; x). By 
making all of these substitutions into the right side of Eq. (19) we find 

3 F 2 (a+e, l+a/2+e/2, e/2; l+a+e/2, a/2+e/2; x) = (l-x)~ e 2 F 1 (a, -e/2; a+e/2+1; x). 

(22) 

Therefore, Eq. (5) holds. 

Proof of Eq. (6). In this case we require an identity for a 3F2 function before 
reducing to 2 F 1 functions. For the right side of Eq. (6) we have 

4 F 3 (a+c?— 1, —a, 2—d—2a, d+2a— 1; d, d+a, 1— d— 2a; y) = 3 F 2 (a+d— 1, —a, d+2a— 1; d, d+a; y) 



+ a(a + d 1) ^ 3F ^ a + dl _ ad + 2a .d+l,d + a + l;y). (23) 
d(a + d) 

We have the contiguous relation [15, 12] 

-a(d + 2a - l)y 3 F 2 (a + d, 1 - a, d + 2a; d + 1, d + a + 1; y) 

= d(d+a) [ 3 F 2 (a + d, —a, d + 2a — 1; d, d + a; y) — 3 F 2 (a + d — 1, —a, d + 2a — 1; d, d + a; y)] 

= d(d + a) [ 2 Fi{—a, d + 2a — 1; d; y) — 3-F 2 (a + d — 1, —a, d + 2a — 1; d, d + a; y)] . 

(24) 

Then to confirm Eq. (4) we use [12] 

c d 
3 F 2 (a, 6, c; d+1, c+ 1; z) = ^—^ 2 Fi(a, b; d+ 1; z) - ^— - 3^2(0, 6, c; d, c+ 1; z), (25) 

and the transformation (18). 

Other 3 F 2 function relations 

Mostly the above reductions depended upon a Gauss or Clausen hypergeometric 
function having a numerator parameter exceeding a denominator parameter by one. 
Here we discuss how this situation can be extended. The conclusion is that if a 
numerator parameter exceeds a denominator parameter by a positive integer k, a 
Clausen 3 F 2 function may be written as a sum of k + 1 2 Fi functions [12]. 

As another example to the previous section, we have 

(a + 2) n n n(n — 1) , , 

= 1 + 2- + -) f, (26) 



(a) n a a(a + 1) ' 

so that 



3 F 2 (6,c,a + 2;d,a;a;) = f; (6)j(c) 



"0 



1 + 2 J_ + J(J-1) 



a a(a + 1) 



x 



be 



6(6 + l)c(c+l) 



2 F 1 (b,c;d;x)+2—x 2 F 1 (b+l,c+l;d+l;x)-[ 

eta a(a + l)a(a + 1) 



x 2 2 Fi(6+2,c+2;rf+2;x). 



(27) 



In general, we have 

(a + k) n (a + n) k 1 



(a)n (a)fc (a) k 



{a)k + (^jn(a + 1)^ + (^jn(n - l)(a + 2) fe _ 2 



+ ...+ n(n-l)--(n-k + l) 



(2f 



a result provable by the method of finite differences. Therefore the general reduction 
is given by 

3 F 2 (b, c,a + k; d, a;x)=j2 rWk** ^ b + + i; d + i; x). (29) 

e=o W WtWt 

Other 3 F 2 function reductions 

The authors of Ref. [6] used identities given in Abramowitz and Stegun [1] to 
partially transform two particular 3 F 2 functions. We complete and then illustrate 
this development. 

We have ([6], p. 1987) 

a4 b r 1 



(74 f L 

3 F 2 (a,b,b + l/2;a + l,2b;z) = — (1 - y)-^ + y)- a dy, (30) 

Z a J(l-z) 1 / 2 



and 



3 F 2 (a, b, b + 1/2; a + 1, 1/2; z) = ~ a f" y 2 "" 1 [(1 + y y 2b + (1 - y)- 2b ]dy. (31) 

z a Jo 

In Eq. (30) or (31), one should take b ^ to avoid a trivial case. In both Eqs. (30) 
and (31) one should take a ^ —1 to have a well defined 3 F 2 function. The form of the 
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integrand factors indicates when certain special cases may arise: for instance when 
a = 1, a = 26, or 6 = 1/2 for Eq. (30) and when a = 1/2 or a = 1/2 - b for Eq. (31). 
Indeed, we give examples of such cases below. It is also obvious that many cases of 
Eqs. (30) and (31) can be written in terms of the incomplete Beta function. 

We may perform each of the integrations in Eqs. (30) and (31) in terms of the 
Gauss hypergeometric function. We do this by introducing a change of variable 
to transform the integral to the interval [0,1], and then apply a standard integral 
representation for 2 F\ (e.g., [1, 2, 7, 12, 17]). 

For Eq. (30), we let v(y) = [y - (1 - z) l / 2 ]/[l - (1 - z) 1 / 2 } and there results 

3 F 2 (a, b, 6+1/2; a+1, 26; z) = 4 b [l+(l-z) 1/2 ]- 2b 2 F 1 [2b - a, 1; a + 1; 1 + -[(1 - z) 1 ' 2 - 

(32) 

We may then use the transformation (e.g., [1, 7, 12, 17]) 2^1(0;, f3; 7; z ) — (1 — 
z)~ a 2 Fi[a,7 — fl;~f;z/(z — 1)]. After some algebraic manipulations, we obtain the 
alternative form 

3 F 2 (a, 6, 6+1/2; a+1, 26; z) = - ^ z)1/2]a 2^1 (26 - a, a; a + 1; ^[1 - (1 - zf' 2 }^ . 

(33) 

For Eq. (31) we simply use v(y) = y/z 1 / 2 and find that 

3 F 2 (a, 6, 6+1/2; a+1, 1/2; z) = - [ 2 F 1 (2b, 2a; 2a + 1; -z 1/2 ) + 2 F ± (2b, 2a; 2a + 1; z 1 ' 2 ) 

(34) 

We have furnished reductions of both Eqs. (30) and (31) for general values of the 
parameters a and 6. 
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Equations (33) and (34) have a multitude of special cases, several of which return 
known results [1, 7]. If a = 1 in Eq. (33) and we apply a known reduction for 2 i*i [7], 
we obtain 

3 F 2 (1, 1 - n/2, 3/2 - n/2; 2, 2 - n; z) = — [l - 2" n (l + (1 - z) 1/2 ) n l . (35) 

nz L J 

In this case, 6 = 1 — n/2 and ra need not be an integer. For the very special case 
that 6 = 1 in Eq. (33), the reduction 2 -f\ (1, 1; 2; — z) = ln(l +z)/z [7] can be applied, 
giving [9] 

3 F 2 (1, 1, 3/2; 2, 2; z) = - fin 2 - ln(l + (1 - z) 1 ' 2 ) 



= -lln 

z 



1 + (1-^)V2 



(36) 



When a = 26 in Eq. (30) or (33) we have 



3 F 2 (26, 6, 6 + 1/2; 26 + 1; 26; z) = 2 Fi(6, 6 + 1/2; 26 + 1; z) = [l - (1 - z) 



vl/2 



26 



)26 



[l + (l-z) 1 /2]26- 

This equation is equivalent to a known result [1]. 
For a = 1/2 in Eq. (34) we have 

3 F 2 (l/2, 6, 6 + 1/2; 3/2, 1/2; z) = 2 F 1 (b, b + 1/2; 3/2; z) 
1 1 



(37) 



(1 + z l/2 ) l-26_ (l _ z l/2 ) l- 



26 



(38) 



2zV2 (i _ 26) 

that is also equivalent to a known result [1]. When 6 = 1 in this equation, we have 
the very special case 2 iq(l, 3/2; 3/2; z) = i-F (l; —;z) = 1/(1 — z). We may inquire 
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as to the case of b — > 1/2 in Eq. (38), and in so doing make use of the expansion in 
powers of 1 — 2b resulting from writing (1 ± z l l 2 ) l ~ 2b = exp[(l — 2b) ln(l ± z 1 ^ 2 )]. The 
result is 



2^(1/2,1; 3/2; z) 



1 , fl + z 1 / 2 ' 
In 



2^/2 I 1 - z l/2 



A/2 



In 



l + z 



1/2 



{i- z y/ 2 



, 0<z<l, (39) 



equivalent to a previous result obtained in a different way [16]. For 5 = a+ l/2in 
Eq. (34) we obtain 

3 F 2 (a, a + 1/2, a + 1; a + 1, 1/2; z) = 2 F 1 (a, a + 1/2; 1/2; z) 

2 F l (2a + 1, 2a; 2a + 1; -z 1/2 ) + 2 ^i(2a + 1, 2a; 2a + 1; z 1/2 )] 



1 F (2a; -; + lFo(2 a; -; -z 1 ' 2 )] = + z 1 / 2 )^ + (1 - z 1 ' 2 ) 



a known result [1]. When b = a — 1/2 in Eqs. (31) and (34) we have 

1 



-2a 



(40) 



3 F 2 (a, a-1/2, a; a+1, 1/2; z) 



Fx (2a - 1, 2a; 2a + 1; -z 1/2 ) + 2 Fx(2a - 1, 2a; 2a + 1, z 1/2 ) 

(41) 

By the standard transformation (18) (e.g., [1, 2, 7, 17]) the hypergeometric functions 
in this equation are given by 2 -Fi(2a— 1, 2a; 2a+l; w) = (l—w) 2 ^~ a ^ 2^(2, 1; 2a+l; w). 
Other cases include b = 1/2 — a in Eq. (31) and b = 1/2 in Eq. (30), and still 
others may be considered. Our examples well illustrate the range of 3F2 reductions 
obtainable from Eqs. (30)-(34). 
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Concluding remarks 

We have given reductions of some particular single variable hypergeometric func- 
tions. These and related series appear in a variety of physical science applications 
and elsewhere, including the power series solution of the Schrodinger equation and 
the summation of perturbation expansions in quantum mechanics. Our direct proofs 
of the reductions should serve to illustrate many of the properties of the functions 
q F q and q +iF q . More familiarity with such series as we have considered may allow 
physicists, quantum chemists, and others to more often use them and to potentially 
increase the computational efficiency of numerical methods, for instance, in many- 
body wavefunction calculations. 

The symbolic implementation of hypergeometric function identities has become a 
more and more active area of research [2, 10, 11]. When a sum can be written as a 
hypergeometric series, many techniques can be brought to bear, including the W-Z 
(Wilf-Zeilberger) method [2, 10]. As part of such methods, a recurrence relation is 
written that shows a "creative telescoping'" [2, 10]. Hence new symbolic computing 
tools [10, 11] are becoming available. 
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